Dynamical zero-temperature phase transitions and cosmic inflation/deflation 
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For a rather general class of scenarios, sweeping through a zero-temperature phase transition by 
means of a time-dependent external parameter entails universal behavior: In the vicinity of the 
critical point, excitations behave as quantum fields in an expanding or contracting universe. The 
resulting effects such as the amplification or suppression of quantum fluctuations (due to horizon 
crossing, freezing, and squeezing) including the induced spectrum can be derived using the curved 
space-time analogy. The observed similarity entices the question of whether cosmic inflation itself 
might perhaps have been such a phase transition. 
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In contrast to thermal phase transitions occurring 
when the strength of the thermal fluctuations equals a 
certain threshold (and so changes the character of the 
stable phase) , zero-temperature phase transitions such as 
quantum phase transitions Q] denote the crossover of dif- 
ferent ground states at a certain (critical) value of some 
external parameter, where quantum fluctuations play a 
dominant role, cf. Fig.^ In both cases, there exists a vast 
amount of literature regarding the equilibrium properties 
in the vicinity of the phase transition, for example in view 
of universal behavior (e.g., scaling laws) near the critical 
point. However, since response times typically diverge in 
the vicinity of the critical point, sweeping through the 
phase transition with a finite velocity, for example, leads 
to a break-down of adiabaticity and thus might gener- 
ate interesting dynamical (non-equilibrium) effects. For 
thermal phase transitions, a prominent example is the 
Kibble-Zurek mechanism, i.e., the generation of topolog- 
ical defects via rapid cooling (quench), which can be ap- 
plied to the phase transitions in the early universe as well 
as in the laboratory 
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FIG. 1: Sketch of level structure near the critical point. 
The energy E of the levels is plotted as a function of some 
external parameter g. At the critical point g — g c , the ground 
state changes from |5 r <) to |5 r >). Typical excitations above 
the ground state |^<) for g < g c are denoted by $ and \. 
Some of the excitations x may remain stable after crossing 
the critical point g = g c , whereas others <E> become unstable 
i.e., lie below | v I / <}. 



The following considerations are devoted to non- 
equilibrium effects in a rather general class of zero- 
temperature transitions, which also expose a remarkable 
analogy between cosmology and laboratory physics, see 
also [3| . Let us consider a quantum system (at zero tem- 
perature) described by the Hamiltonian H depending on 
some external parameter g. At a certain critical value of 
this parameter g c , the system is supposed to undergo a 
phase transition, i.e., the ground state \ 1 $ < (g)) of H(g) 
for g < g c is different from the ground state \^>{g)) 
of H(g) for g > g c . For example, \9 < (g)) and \^f > (g)) 
could have different global/ topological properties (such 
as magnetization) in the thermodynamic limit. 

Before the phase transition g < g c , the phase space can 
be explored in terms of quasi-particle excitations above 
the ground state |\E , <(g)). If we make the assumption 
that the Hamiltonian H is analytic in the external param- 
eter g, we may follow the behavior of these excitations 
through the critical point g c and extrapolate them to the 
region g > g c . For g > g Cl the state \^<(g)) is no longer 
the ground state and hence some of the excitations must 
become unstable after the critical point g > g c , since now 
\^>{g)) has a lower energy than |\E' < (5)) , cf. Fig. ^ In 
the following, we shall focus on these quasi-particle exci- 
tation modes, whose energy is positive for g < g c and be- 
comes negative after crossing the critical point g > g c . In 
order to describe these excitations quantitatively, a few 
additional assumptions/approximations are necessary: 
Linearity: We assume that the crossover from stability 
to instability of these modes can be described using lin- 
ear stability analysis, i.e., that the associated quantum 
fluctuations are small enough. Consequently, the onset 
of instability occurs if the dispersion relation uj 2 (k) dives 
below the fc-axis. 

Vanishing gap: Since we are mainly interested in low- 
energy and long-wavelength excitations (which will turn 
out to yield universal behavior largely independent of the 
microscopic structure) , we shall assume that the modes $ 
are gap-less such as Goldstone modes, i.e., uj 2 (k = 0) = 
for all g. The analogy to quantum fields in an expand- 
ing/contracting universe applies to the general case with 
a gap as well, but the main points of interest (analogue 
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of cosmic horizon) can be studied for gap-less modes. 
Analyticity: The dispersion relation u> 2 (k,g) is sup- 
posed to be an analytic function of k. Hence a Taylor 
expansion starts with the term [with c 2 {g < g c ) > 0] 

u; 2 (k,g) = k 2 c 2 (g)+0(e). (1) 

Independence: As the final ingredient, we assume that 
the unstable modes are independent of each other such 
that it suffices to consider one scalar mode <£>. 

Based on these assumptions, we may construct the 
low-energy effective theory describing the mode $ in the 
vicinity of the critical point g c according to the disper- 
sion relation QJ. As demonstrated in Ref. Q, given the 
above conditions, the most general low-energy effective 
action can be cast into the form 

£ eff = \ g£ (d, $) (d u <f>) , (2) 

where \f—g c s g^s denotes a matrix depending on the sys- 
tem under consideration as well as the external parameter 
and a sum convention over [a, v = . . . 3 is employed. 

Hence the excitations $ behave as minimally coupled 
and mass-less scalar quantum fields in curved space-times 
whose geometry is determined by the effective metric 
g^g', which is basically the underlyingidea of the ana- 
logue gravity concept, see, e.g., [a, la 0- If we assume 
the quantum system under consideration to be effectively 
homogeneous and isotropic at large wavelengths A, the 
above action simplifies to 

£ cff = i(i$ 2 -/3(V$) 2 ) . (3) 

The external parameters a and (3 may depend on g and 
hence on time and must be non-negative before the phase 
transition. The convenience of the choice of 1/a in the 
first term becomes apparent after constructing the asso- 
ciated effective Hamiltonian 

Weff = ^(an 2 +/3(V$) 2 ) ■ (4) 

In the homogeneous and isotropic case, the effective met- 
ric reads (for g < g c ) 

ds 2 cS = JafP dt 2 - yfWfadr 2 . (5) 

Since the mode <& becomes unstable for g > g c , at least 
one of the two parameters has to change its sign at the 
critical point. There are basically three possibilities: 
A: q | while (3 remains finite. (Note that 1/a J, is 
not possible since the propagation speed would diverge.) 
According to Eq. (JSJ , the effective metric corresponds to 
an expanding universe in this case. 

B: j3 I while a remains finite (similarly, 1/(3 [ is not 
possible), which corresponds to a contracting universe. 
C: Both a and (3 become singular. 



Having established the analogy to cosmology, we may 
now apply the tools and concepts known from curved 
space-times 8J. The space-time described by the above 
metric contains a (particle) horizon if the maximum dis- 
tance (measured in co-moving coordinates) which can be 
traveled starting at the time t- ln 

tout tout 

Ar=y dtc(t)= J dt y/a(t)/3(t) , (6) 

tin tin 

is finite, i.e., points beyond the horizon size Ar can never 
be reached. Obviously, since c = yja(t)/3(t) is bounded 
from above, a horizon always exists if the critical point 
is reached in a finite laboratory time t out (which might 
still correspond to an infinite proper time in cosmology), 
but also if a and/or (3 decrease fast enough for £ out | 00 ; 
for example a a l/t n with n > 2, cf. 

As one can infer from the above equation, the horizon 
size (in terms of the laboratory coordinates t, r) always 
decreases as ti n increases. Hence, all $-modes with a 
given (effective) wavelength A cross the horizon at some 
time (when A exceeds the horizon size Ar). After this 
horizon crossing, the modes cannot oscillate anymore 
(freezing) due to loss of causality across the horizon and 
their quantum state gets squeezed. This mechanism is 
completely analogous to the amplification of quantum 
vacuum fluctuations of the inflaton field in our present 
standard model of cosmology - which are supposed to be 
the seeds for structure formation. 

In order cast this analogy into a more quantitative 
form, let us consider the equations of motion for $ 

In the homogeneous and isotropic case Va = V/3 = 0, 
there is a duality between the field $ and its canonical 
momentum n = $/a since n obeys the same the equa- 
tion of motion as <& with a and (3 interchanged 

(l^)|-^ V2 ) n ^°- (8) 

This duality is analogous to that between the electric 
and magnetic field in the absence of macroscopic sources. 
However, even though the equations of motion in cases 
A and B are related by this duality, the behavior of the 
quantum fluctuations (e.g., their spectrum) is different. 

Assuming general power-law time-dependence near the 
critical point a cx t a and (3 oc t b , we may remem- 
ber the coordinate invariance of the effective geometry 
and introduce another time coordinate r proportional to 
i 2 / ( 2 + a + b ) j for which the wave equation simplifies to 

(£-^!-*)-o. <*) 
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provided that v is given by 

1 + a 



(10) 



Hence the solutions for the spatial Fourier modes can be 
expressed in terms of Hankel functions |lfll | 



d>±(r)=r^±(fcr) 



(11) 



Sufficiently far away from the critical point |r| f oo, 
the modes oscillate 3?iT(t) ~ t" e ±tkT /\/~kr and the so- 
lutions $^ are basically the functions in the expansion 
into creation and annihilation operators corresponding to 
the adiabatic vacuum state, cf. Q. When approaching 
the phase transition, however, adiabaticity breaks down 
and the modes cross the horizon and freeze [lfj 



*fc(M I o) ~ r 



(12) 



I.e., the spectrum of the two-point correlation function 
($(r)$(r')) of the frozen quantum fluctuations is deter- 
mined by the parameter v > in Eq. i|l(J|) . 

Let us examine a few examples: The trivial case 
a = b = of course reproduces the undisturbed spectrum 
v = 1/2. Neglecting the back-reaction of the quantum 
fluctuations onto the the dynamics of the external pa- 
rameter g, its time-evolution g(t) should not experience 
anything special at the critical point g = g c and hence 
the most natural choice for its dynamics is a constant 
velocity g — g c oc t. If we assume the effective Hamilto- 
nian to be an analytic function of g, this implies a = 1 
and/or 6=1. (Otherwise a and 6 are the character- 



istic exponents \g — g c \ a and \g 



occurring in H .) 



In case A (expanding universe) with a — 1 and 6 = 0, 
we obtain v = 2/3, i.e., quantum fluctuations with large 
wavelengths are amplified, case B (contracting universe) 
with a = and 6=1 yields v= 1/3, i.e., quantum 
fluctuations with large wavelengths are suppressed. Fi- 
nally, a = b = 1 (case C) leads to an undistorted spec- 
trum v = 1/2, which is not surprising once one realizes 
that the effective metric in Eq. (JSJ is exactly flat in terms 
of the new time-coordinate ds 2 s = dr 2 — dr 2 . 

Now we are in the position to apply the above method 
to some concrete physical systems: As a first example, we 
consider atomic Bose-Einstein condensates. In the dilute- 
gas limit, the quantum phase and density fluctuations 
are small and can be treated as linear perturbations. For 
wavelengths far above the healing length, the effective 
action of the phase fluctuations <& reads (H = 1) 



1 fl 



Lh) 



c eS = ~ - & - ^ (V$) 

g m 



(13) 



where go is the background density of the condensate, m 
the mass of the atoms, and g the time-dependent coupling 
strength representing the inter-particle repulsion g > or 
attraction g < 0. Obviously, a homogeneous condensate 



becomes unstable for attractive interactions (i.e., g c = 0) 
and this scenario g J. corresponds to case A, i.e., an 
expanding universe. Hence sweeping through the phase 
transition at g c = by means of a time-dependent exter- 
nal magnetic field (Feshbach resonance) with a finite ve- 
locity generates a fc~ 4 / 3 -spectrum of the two-point phase- 
phase correlation function. However, it might be difficult 
to measure the phase after the phase transition (collapse 
of condensate due to molecule formation), whereas the 
frozen density fluctuations with large wavelengths should 
easily be measurable. The density fluctuations are just 
the canonically conjugated momentum field n and can be 
calculated analogously using the duality in Eqs. J7J and 
© , but with an additional factor of k in Eq. ifTTl) . Hence 
the spectrum of the density-density correlation function 
behaves as fc +4 / 3 . Note that this behavior is consistent 
with the amplification/suppression of quantum fluctua- 
tions by squeezing which maintains the minimal Heisen- 
berg uncertainty of the ground state, i.e., Aqk Apk = h/2. 

As a second example, let us study a simple 1+1 di- 
mensional model of the electromagnetic field coupled to 
a linear medium via the magnetic component 



Ccff = 2 ( 



E 2 - B 2 



* 2 - n 2 



* 2 + gB*) , (14) 



with the electromagnetic field being governed by the po- 
tential A via E = d t A and B = d x A. The field * de- 
scribes the (linearized and localized) dynamics of the 
medium with the plasma frequency f2 and g denotes the 
coupling (magnetic dipole moment). Averaging over the 
degrees of freedom ^ of the medium, the low-energy ef- 
fective theory for the macroscopic electromagnetic field 
yields the permeability = 1 — g 2 /£l 2 (which corre- 
sponds to inserting the adiabatic solution ^ « gB/fl 2 
back into the action). Hence there is a phase transition 
at the critical value of the coupling g c — Q after which 
the medium develops a spontaneous magnetization and 
the linearized description above breaks down. Identifying 
A = $, this scenario corresponds to case B. Hence the 
frozen two-point spectra behave as k~ 2 / 3 for A, and thus 
£,4/3 f or q ^ anc j an( j g na Uy £2/3 f or £ _ jj ( a g a i n re- 
specting AqkApk — H/2) if we sweep through the critical 
point with a finite velocity. Note that, for the source- 
free macroscopic electromagnetic field, one can also in- 
troduce the dual potential A via H = B/fi = 9 t A and 
D = E = d x A, which explicitly incorporates the duality 
in Eqs. (J7J) and (JSJ, i.e., case B — > case A, and conse- 
quently has a fc~ 4 / 3 -spectrum. 

Finally, as a third example, we consider a very simple 
version of the Heisenberg model 



H 



-g(t) 



<ij> 



(15) 



with (Ti being the Pauli spin-1/2 operators for each lat- 
tice site i and < ij > denoting the sum over nearest and 
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next-to-nearest neighbors, for example. The ferromag- 
netic state |^<) = |. . . tttt ■ • • ) is the ground state for 
g > and breaks the 0(3) invariance of the Hamiltonian; 
thus the spin-waves (magnons) are gap-less (Goldstonc 
theorem). Obviously g = is the critical point here, 
after which the energy of the magnons becomes nega- 
tive. In view of the global factor g(t), this scenario is 
an example for case C with a = /3 = g, i.e., a = b. The 
observation that the spectrum of the fluctuations is not 
disturbed by the dynamics of g(t) can also be explained 
by the fact that the ferromagnetic state |... TTTT - -- ) 
is an exact eigenstate of the Hamiltonian, i.e., even the 
decay from |\&<) to |\&>) is impossible with the exact 
Hamiltonian in Eq. (|15JI and requires some disturbances. 

In summary, the analogy between the excitations <I> 
and quantum fields in an expanding/contracting universe 
allows us to apply universal geometrical concepts such as 
horizons to a rather general class of quantum systems 
approaching the critical point. Near the phase tran- 
sition, adiabaticity breaks down (since the energy gap 
vanishes) and the system does not stay in its (instanta- 
neous/adiabatic) ground state in general. The spectrum 
of the two-point correlation function of the quantum fluc- 
tuations frozen out at the phase transition (which can be 
the seeds for pattern formation etc.) can be derived quite 
independent of the microscopic details of the considered 
system and is basically determined by the characteristic 
exponents a and b (universal behavior). The strength of 
the frozen fluctuations, however, and their dynamics af- 
ter crossing the critical point (nonlinear instabilities etc.) 
depend on the explicit microscopic structure (e.g., the di- 
luteness parameter in Bose-Einstein condensates). 

As an outlook, one might compare phase transitions 
(such as the examples considered above) to "real" cos- 
mic inflation which is part of the present standard model 
of cosmology. Interestingly, phase transitions reproduce 
many features of inflation: The decay of |^<) down 
to |^>) at g > g c (breakdown of adiabaticity) releases 
energy (analogous to reheating). Phase transitions dis- 
play universal behavior (no fine-tuning) in the sense that 
initial small-scale deviations from | V I'<) are not impor- 
tant after the transition g > g c . Similarly, the large- 
scale homogeneity may be explained naturally. How- 
ever, even though quantum fluctuations generate small 
inhomogeneities in both scenarios (phase transitions and 
"real" cosmic inflation), none of the examples consid- 
ered above reproduces the correct scale-invariant 1/fc 3 - 
spectrum. However, that is not surprising as the ex- 
amples considered above break many symmetries we ob- 
serve in the real universe, e.g., they possess a preferred 
frame and do not respect the principle of equivalence etc. 
If we demand that the effective action (at least at low 
energies) does not single out a locally preferred frame 
(e.g., the two-point function (&(x)$(x!)) depends on the 
Ricci scalar etc.) there are only two possibilities: firstly, 
the effectively flat space-time with a = b leading to an 



undisturbed 1/fc-spectrum, and, secondly, the remaining 
nontrivial combination a + 36 = —4 exactly corresponds 
to the de Sitter metric and thus reproduces the scale- 
invariant l/fc 3 -spectrum, cf. Eq. IjlUfl and 9J. The sec- 
ond version allows for more additional symmetries: We 
may assume a constant velocity of propagation for the 
3>-mode 

f , 1 <j> 2 - (V$) 2 
A=jdtd 3 r- ^J-, (16) 

and the resulting action turns out to be scale- invariant: 
A[t — > Xt, r — * Xr] = A[t, r]. (There are also other moti- 
vations for the above form of the action such as the prin- 
ciple of equivalence, which shall not be discussed here.) 
One would expect the dynamics of this action, which has 
been motivated by demanding the above symmetries, to 
be generated by the back-reaction, which has been omit- 
ted so far and which respects these symmetries. These 
interesting findings entice the question/speculation of 
whether cosmic inflation itself might perhaps have been 
such a phase transition. 
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